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ABSTRACT 

The arm structure of the Milky Way remains somewhat of an unknown, with observa¬ 
tional studies hindered by our location within the Galactic disc. In the work presented 
here we use smoothed particle hydrodynamics (SPH) and radiative transfer to create 
synthetic longitude-velocity observations. Our aim is to reverse-engineer a top down 
map of the Galaxy by comparing synthetic longitude-velocity maps to those observed. 
We set up a system of A-body particles to represent the disc and bulge, allowing 
for dynamic creation of spiral features. Interstellar gas, and the molecular content, 
is evolved alongside the stellar system. A 3D-radiative transfer code is then used to 
compare the models to observational data. The resulting models display arm features 
that are a good reproduction of many of the observed emission structures of the Milky 
Way. These arms however are dynamic and transient, allowing for a wide range of 
morphologies not possible with standard density wave theory. The best fitting models 
are a much better match than previous work using fixed potentials. They favour a 
4-armed model with a pitch angle of approximately 20°, though with a pattern speed 
that decreases with increasing Galactic radius. Inner bars are lacking however, which 
appear required to fully reproduce the central molecular zone. 

Key words: hydrodynamics, radiative transfer, ISM: structure. Galaxy: structure, 
kinematics and dynamics, galaxies: spiral 


1 INTRODUCTION 

The question of the structure of the Milky Way is still an 
open one. A top-down view of our Galaxy has been sought 
after for decades, and yet due to our location within the 
disc itself such a map remains elusive. The general consensus 
favours a multi-arm spiral pattern with some kind of inner 
bar structure, though the exact morphology and kinematics 
of these features are less well constrained (Sewilo et al. 2004| 


Vallee|2005| |Benjamin|2008| [Francis fc Anderson|2012 |Pet-| 

titt et al.|2014 ). For example, spiral models with 2 or 4 arms 

have been proposed to best represent Galactic morphology 
(Elmegreen 1985 Liszt 1985| Georgelin & Georgelin 1976 


Taylor & Cordes 1993| | Weaver || 1970 1 ) though there is evi¬ 
dence that the gaseous and stellar components trace out dif¬ 
ferent arm numbers, with gas favouring a 4-armed structure 
and stars a 2-armed one ( Drimmel|2Q00| Drimmel & Spergel 


2001 Benjamin 2008 Steiman-Cameron et al. 2010). Fur¬ 


thermore, the pitch angle of these arms also varies dramat¬ 
ically between studies, with a preference to tighter wound 
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arms in the inner disc compared to the outer, with values 
ranging from 5° to 20° ( Levine et al.|200^ Hou et al.|2009 


Efremov]|2011 Francis Anderson||2012'). 

The status of the Galactic bar further complicates 
the already confusing picture. Values for the orientation, 
bar length, and rotation speed vary between studies, es¬ 
pecially the rotation speed for whi ch values range from 
30kms“^kpc“^ to 70kms“^kpc“^ (Gerhard 2011). The 


suggestion of a second narrow bar in addition to the clas¬ 


sical 3kpc long COBE DIRBE bar ( Blitz k, Spergel||l991 

Weiland et al.[19^ Binney et al.|1997 Gerhard 2002| only 


confuses matters further ( Hammersley et al.|2000 Benjamin 

et al.||2005 ). 

Determining the distances to sources to map out Galac¬ 
tic structure has many difhculties and uncertainties. Using 
longitude-velocity maps {l-v) can circumvent this issue by 
displaying spiral arm features without need of calculating 
distances. Such maps exist in a number of different sources 


( Kalberla et al.|2005 Strasser et al.|2007 McClure-Griffiths 


et al. 

2012 

Caswell & Haynes] 1987|[Cohen et al.|1980|[Dame 

et al. 

2001 

Jackson et al.|2006| [Pineda et al.|2013[). The CO 


l-v emission map of Dame et al. (2001) is shown in Fig- 
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Figure 1. Longitude-velocity map of brightness temperature of 
the CO (J=0-1) transition ( |Dame et al.||200T] >, with major arm 
features labelled. We integrate the CO emission over latitude in 
order to show weaker features. Ql-4 indicates the position of 
Galactic quadrants. 


ure[^ where molecular gas can be seen to trace out the spiral 
structure as regions of increased emission strength. 

Top-down maps of the Galaxy can be recreated by 
reverse-engineering models to match the features seen in 


1999 

|G6mez & Gox|[2004 

IRodriguez-Fernandez & Gombes 

2008 

|Baba et aL| 

20101 ii 

“Choperskov et al. 2013[ Sormani 


titt et al. (2014), hereafter referred to as Paper 1, we used 


fixed analytical gravitational potentials for the arm and bar 
components in simulations of ISM gas to constrain Galac¬ 
tic structure. Synthetic longitude-velocity maps were con¬ 
structed from simulations to compare directly with the ob¬ 
served GO emission. A wide array of potential morphologies 
were tested, including numerous pitch angles, pattern speeds 
and arm numbers. We found arm and bar features were nec¬ 
essary to match the observations, with arm pattern speeds 
of 20km s“^ kpc“^ and bar pattern speeds in the range of 50- 
60km s“^ kpc“^. The 4-armed models could reproduce many 
of the features seen in observations but produced too much 
excess emission. The 2-armed model provided too few fea¬ 
tures compared to observations. The symmetric nature of 
such potentials, used to represent the spiral density waves 
of Lin & Shu (1964), made the exact placement of arms 


in l-v space problematic. For example, no 4-armed model 
could satisfactorily reproduce the Garina and Perseus arms 
simultaneously. 

In this paper we take a different approach and model 
the stellar distribution of the Milky Way with a set of dis¬ 
crete Wbody particles rather than a continuous potential 
as in Paper 1. This method has been used to simulate both 
isolated stellar systems ( Sellwood Garlberg||1984[ Sparke 


&: Sellwood|1987l|Sellwood Sparke|1988||Shen et ai.|2010[ 

Grand et al.||2Q13f D’Onghia et al.||2Q13 ) and the simulta- 


neons evolution of a gas disc ( 

Garlberg & Freedman 


Elmegreen & Thomasson 

1993| 

Glarke & Gittins|20061 

Baba 

et al. 

2009 

Grand et al. 

12012 

[Struck et aL||2011| Renaud 

et al. 

2013 

Athanassoula et aL||20131 Roca-Fabrega et al. 


shown that both bar and spiral features are surprisingly easy 
to reproduce, though not ones that necessarily agree with 
theory. For instance, spiral arms do not appear to be steady 
spiral density waves such as those suggested by |Lin Shu| 
(1964), and bars are seemingly overabundant, and difficult 


to reproduce with accompanying spiral structures. Instead, 
properties such as number of arms and pattern speeds are in 
better agreement with swing amplified instabilities ( [Toomre 


1981 


Grand et aL||2Q13 see Section 2.3). 


This paper is organised as follows. In Section we out¬ 
line the numerical methods used in this paper and the setup 
procedure for the various Galactic components. In Section 
we present the results of our simulations. This is divided into 
the hydrodynamical simulations (|3.1|and 3.2) and l-v maps 


constructed using a simple approximation (3.3) and a full ra¬ 
diative treatment (3.4). In Sectionj^we discuss the successes 
and limitations of this study, specifically with comparison to 
the analytical potentials of Paper 1, and we conclude in Sec¬ 
tion [b] 


2 NUMERICAL SIMULATIONS 

We use the SPH code SPhNG based on the original version of 
Benz et al. (1990), but substantially modified as described 
Bate et al. (1995) and Price & Monaghan (2007), and 


parallelised using both OpenMP and MPI. This is the same 
code as was previously used to study the spiral structure of 


M51 by Dobbs et al. (2010). 


2.1 Hydrodynamical simulations 

The stellar distribution is represented by a collection of A- 
body particles constituting the disc and bulge components, 
and in some cases the halo. Star particles have adaptive 
gravitational softening lengths, which are calculated in a 
similar manner to the smoothing lengths for SPH particles 
( Price Monaghan||2007 ). 

The gas disc itself is devoid of self-gravity, where we 
assume that the majority of the gravitational force stems 
from the stellar material. Gas self-gravity may dominate the 
smaller scales, but we are interested predominately in the 
global gaseous structure. As with Paper 1, stellar feedback 
and magnetic fields are not included. The lack of feedback 
will cause gas to collapse into the plane rapidly, though the 
inherent dispersion of the stars keeps the gas disc thicker 
than when using an analytic potential. 

A basic chemical network is also included, identical to 
that of Paper 1, in which the abundances of HI, H 2 and GO 
are evolved for each SPH gas particle. These abundances are 
then used as inputs to the construction of l-v maps. The gas 


is also subject to various heating and cooling effects (Glover 
fc Mac Low|2007 ). 


2.2 Initial conditions and Galactic potentials 

In our calculations we seed the initial positions and veloci¬ 
ties of the star particles from a set of density profiles tailored 
to match the rotation curve of the Milky Way. The Galaxy 
is decomposed into halo, disc and bulge components, the 
density profiles for which are used to set the initial posi¬ 
tions of the star particles by the tabulation and sampling 
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of mass distribution functions. We initialise the velocities 
using the procedures from Hernquist (1993), also used by 


Baba et ah (2009), which primarily requires integrating the 


separate moments of the collisonless Boltzmann equation. 
Velocities of the disc are circular with some dispersion, while 
the velocities of the halo and bulge are set in random orien¬ 
tations. 

Our various calculations are described in Table in Our 
four fiducial calculations are named Ba, Bb, Be and Bd, 
and use a live stellar bulge and disc in a static halo poten¬ 
tial to reduce computational expense (the a-d denoting the 
decrease in disc to halo mass ratio). The Db and Dc models 
do not contain a bulge component and as such are highly 
susceptible to bar formation. The Hb model is the same as 
Bb but with the replacement of a static halo with a live one. 

We choose a NFW profile ( Navarro et al.|19'% ) to repre¬ 
sent the dark matter distribution, which has a density profile 
of the form 


Ph{r) = 


_ Ph,0 _ 

r/rh{l + r/rhY 


(1) 


where ru defines the halo viral radius and ph,o is the density 
scale factor and is given by 


Ph,Q — 


Mh _ 

ln(l + Cnfw) 


^NFW _ 

+ Cnfw ! (1 + Cnfw) 


( 2 ) 


We use a fixed concentration and halo scale length for all cal¬ 
culations of Cnfw = 5 and r 2 oo = 122 kpc. When using a 
live halo it is necessary to truncate the halo at some distance, 
where we truncate the halo distribution using a double ex¬ 
ponential decay with a scale length of 50 kpc. The stellar 


disc is seeded from an exponential density profile (Binney & 
Tremain^|1987 ) 


Pd{r, z) 


Md 

4irR‘^Zd 


exp (-r/rd)sech^(z/zd) 


(3) 


where ra and Zd are radial and vertical scale length, which 
we set to Td = 3.0 kpc and Zd = 0.3 kpc. The stellar bulge is 
represented by a spherical Plummer profile of the form 


Pb{r) 


3Mb rl 

(j-2 _)_ ,.2^5/2 


(4) 


which is also truncated similarly to the halo at a distance of 
4kpc. The bulge scale length is rb = 0.35 kpc and the mass 
is Mb = 1.05 X 10^“ Mq. 

We set gas on similar orbits and positions to the star 
particles, only with much lower masses in accordance with a 
gas disc of the same mass used in Paper 1 for consistency (8x 
10® M©). Gas was not distributed according to the observed 
surface density profile however, but rather in an exponential 
disc. The gas profile is driven by the stellar disc and settles 
after a few lO’s of Myrs. 

The resulting resolutions adopted for each set of com¬ 
putations is approximately 1 million disc star particles, and 
0.1 million bulge particles. A resolution of 3 million gas par¬ 
ticles was used, a compromise between the 5 million used in 
Paper 1 and 1 million. The latter was seen to be too weak 
to produce sufficient emission due to the poorer resolution 
of the high density and CO rich gas. 


2.3 Stability and spiral modes 

The morphology of the stellar disc can be quantified using 
some simple parameters. The Toomre parameter in the stars, 
Qs ( Toomre||1964 ), characterises the stability of the disc to 
local collapse and is given by 


Qs = 


I^O-R 

3.36GEo 


(5) 


where k is the epicycle frequency, an the radial velocity 
dispersion and Eo the disc surface density. Values of Qs < 1 
imply the disc is gravitationally unstable. Physically, Qs can 
be thought of as the balance between the disc pressure forces 
(driven by velocity dispersion in stars or thermal dispersion 
in gas) and gravitational attraction of the disc. The initial 
value of Qs must be dehned for the disc to calculate the 
initial velocity dispersion. We use a value of Qs = 1 to ensure 
the disc is borderline stable to arm formation. 

The dominant swing amplified mode of the stellar disc, 
denoted m ( Toomre|1981 ), can be calculated by 

~ 27rGEoX ~ 47rGEo 

where 1 < V < 2 generates spiral features, and V = 2 


is a nominally adopted value (Fuji! et al. 2011 Dobbs V 
Baba|[20T4 ). Equation]^ shows that the number of arms is 


a strong function of the disc mass, with systems with high 
disc-to-halo mass ratios forming only a few strong spiral 
arms, whereas low mass discs form numerous but weaker 


arms (Carlberg V Freedman 1985). The dominant arm mode 


is also directly coupled to Qs, and so seeding a disc with a 
very low Qs will lower the expected m. Equation E only 
predicts the dominant mode, and any single simulation will 
have other spiral modes of weaker or comparable presence. 
The value of m also increases with radius, and so a simula¬ 
tion should show a greater number of arms in the outer disc 
compared to the inner disc. 

As the disc mass drives the dominant spiral mode (Eq. 
we use four separate disc to halo mass ratios, the rotation 
curves for each are shown in Eigurej^ with the respective 
parameters for the “light”, “normal”, “heavy” and “huge” 
configurations (ordered by increasing disc mass) given in 
Table[^ We keep the bulge initial conditions the same for 
all calculations and use different disc and halo masses to 
reproduce the observed rotation curve. The values for the 
“normal” setup are based on those from Baba et al. (2009). 


The different disc masses also allow for different swing 
amplified spiral modes. While we do not have direct control 
over the spiral structures formed in these calculations, we 
attempt to drive a range of spiral modes induced by swing 
amplification. Table[^ gives the predicted dominant spiral 
mode at three different radii corresponding to the inner, 
mid and outer disc (5, 7.5, and lOkpc). A wide range of 
arm numbers is predicted across all the models, with the Bd 
model aimed at being effectively flocculent. 


2.4 Emission maps 

We use the same method as outlined in Paper 1 to construct 
synthetic emission maps in CO using the radiative transfer 
code TORUS ( Harries|[200Q ). TORUS creates cubes of bright¬ 
ness temperature as a function of latitude, longitude and 
line-of-sight velocity which are then integrated over latitude 
































4 A. R. Pettitt, C. L. Dobbs, D. M. Acreman & M. R. Bate 


Calculation 

Md[lOiOM0] 

M;,[1O1OM0] 

MtllOiOMe] 

m(Rin) 

'^(Rmid) 

'ni(Rout) 

Notes 

Ba 

5.3 

44 

1.05 

2.0 

2.8 

4.0 

“Huge” disc mass 

Bb 

4.1 

63 

1.05 

2.5 

3.6 

5.7 

“Heavy” disc mass 

Be 

3.2 

83 

1.05 

3.1 

4.8 

8.2 

“Normal” disc mass 

Bd 

2.5 

101 

1.05 

4.0 

6.5 

11.5 

“Light” disc mass 

Db 

4.1 

63 

- 

2.2 

3.4 

5.5 

No bulge 

Dc 

3.2 

83 

- 

2.8 

4.6 

8.0 

No bulge 

Hb 

4.1 

63 

1.05 

2.3 

3.3 

6.4 

Live halo 


Table 1. Description of different live-disc models. Refer to Figure for the resulting rotation curves for the main models. If the halo is 
not listed as live then it is represented by an analytic potential. The gas component of the disc is identical in all models with a total 
mass of 8 X 10^ Mq. The swing amplified spiral mode from Equation[^is given at three radii. These are Rin = 5kpc, -Rmid = 7.5kpc and 
Rout = lOkpc corresponding to the inner, mid and outer disc. 



R [kpc] R [kpc] 


Figure 2. Axisymmetric rotation curves are shown for the calculations presented in this study. Disc, bulge and halo components are 
shown as the red, green and blue dashed lines respectively. Four different models are presented with labels referring to the disc mass; 
“light”, “normal”, “heavy” and “huge”. The halo mass is adjusted in each case to ensure a rotation curve that is in agreement with 
observations (blue points, from |Sofne|2012| . 


and q uant it ively compared to the CO map of |Dame et al.| 
(2001). A turbulent velocity of 4kms“^ is added to the line 
width to better match the observed features. 

We must also take into account the uncertainty in the 
observer’s position and velocity by varying the parameters 
Cbs, Robs, Vohs (azimuthal and radial position and circu¬ 
lar velocity of the observer in the disc). As the morphology 
is highly time-dependant we must also test multiple time- 
stamps for each simulation. We choose to do so in the range 
of 200-320Myrs of evolution of the gas disc. This is enough 
time to allow the majority of the molecular gas to form 
(found to be t >150Myr) and for the disc to settle into a 
spiral pattern that will persist for up to a Gyr. Note that 
we must search the full range of lohs values, as we cannot 
assume the disc is symmetric as in the grand design case of 
Paper 1. 

A simple method of creating l-v maps to narrow down 
the parameter space is used, which is described in detail 


in Paper 1. Maps in l-v space are made by projecting the 
gas from x-y to l-v space using /obs, Robs and Wbs, and 
an analogue for CO emission (/i,synth) is calculated using 
the distance from the observer to the gas (di) and the CO 
abundance {xi,co) of the SPH particle. We have however 
slightly changed our equation for the synthetic emission map 
approximation by the addition of an e term 

-^i,synth OC Xi,Co/ T • (7) 

Because of the large amount of interarm material there are 
cases in our fitting where the 1/d‘f approximation causes ex¬ 
tremely large intensity due to particles very close to the ob¬ 
server location. This rarely happened in Paper 1, as molec¬ 
ular gas tended to be exclusively placed along the potential 
minima. To rectify this we include a limiting factor, e, in the 
denominator to limit emission from material near the Solar 
position. To find appropriate values of e we constrain a sim- 
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pie map made from Eq. |7]to a map created from the TORUS 
radiative transfer code, which gives a best fitting value of 
e = 1.0. We maintain k — 2 from Paper 1, finding it to be 
equally good a fit as previously. 

At each main time-step after lOOMyrs we created a 
number of simple l-v maps using Equation and differ¬ 
ent values for the azimuthal position (/obs), radial position 
(^obs) and circular velocity (Wbs) of the observer. A mean 
absolute error (MAE) statistic is used to quantify the differ¬ 
ence between the maps and the data in Eigure^ This is a 
sum over the pixels in a map of the absolute difference be¬ 
tween brightness temperature of the synthetic and observed 
CO maps (see Paper 1 for a more detailed description of 
this process). The fit for several neighbouring time-steps is 
compared, providing a best-fitting time, t?,/, for each model. 
The radiative transfer maps built by TORUS are then used 
to quant it ively compare different models and the CO obser¬ 
vations. 


3 RESULTS 

We devote the majority of this section to discussion of our 
fiducial set of models. The SPH simulations and top-down 
maps are presented first, and both kinds of l-v maps after. 


3.1 Simulations and morphology of fiducial 
models 

The main Galactic disc configuration investigated is the Ba- 
Bd group of models, which contain a live disc and bulge, but 
maintain a static dark matter halo. We show the Ba, Bb, Be 
and Bd models at three separate intervals in Eiguresj^ 
and Irrespectively. Stars are shown in upper panels, and gas 
in the lower. The third time-frame of the Ba model is slightly 
earlier than the other models due to the computation taking 
significantly longer than the other^ We did not continue 
this calculation further as maps made from existing time- 
steps and lower resolution calculations that had evolved for 
longer showed the Ba model was a poorer fit to Galactic l-v 
data in general. Eor each model we attempt to quantify the 
arm number (N), pitch angle (a) and pattern speed (f^p), 
the methods for which are described in the Appendix. 

The gas clearly traces the arm features in the stars in 
Eigures[^|^[^and[^ with the higher density gas in the disc 
being the location of the majority of the molecular content. 
In the highest mass disc (Ba) there is very little inter-arm 
structure, with large regions devoid of gas (Eig. [^. In gen¬ 
eral the arm/inter-arm contrast is much stronger in the gas 
than in the stars. The arms in the Ba model appear non- 
logarithmic in many places, especially so in the final panel, 
with various kinks, or knees, forming along the spiral arms. 
These irregular structures also exist in other models, but are 
strongest in Ba. The spiral structure is much more regular 
in the moderate mass discs; Bb and Be. 

The individual Eourier power spectra, fit to pitch an¬ 
gle and pattern speeds for each model are also given in the 


^ Run time appeared to scale with stellar disc mass due to time- 
step limiting based on the accelerations of the gas particles mov¬ 
ing from very low to high density. 


Appendix, with the best-fitting values presented in Table 
We give N and Qp at three different radii, an inner, mid 
and outer region of the disc. The pitch angle is only quoted 
at a single radius as many R-0 pairs are required for the 
determination. The values in Table [2] are calculated at the 
time where the l-v map best-matches the Dame CO map 
(see Sec. 3.3). The arm numbers show a clear correlation 
between dominant arm mode (in bold) and disc mass. The 
agreement between the dominant value of N{R) and the 
inner-mid m{R) values (Table[^ is good, but the direct com¬ 
parison at specific radii is poorer. We infer that Equation 
is a good rough estimate of the expected dominant arm 
number in a simulation, but cannot be expected to exactly 
reproduce the dominant mode at individual radii. 

In all of the live disc-bulge models the dominant arm 
number is seen to increase with radius, similar to the ra¬ 
dial dependance predicted by Equation In Bb for exam¬ 
ple there is a strong inner 2-armed structure, and a strong 
4-armed component in the outer disc (see the Eourier spec¬ 
trum of Eig. A2 and the top down maps in Eig.[^. The Milky 
Way is thought by some to have a higher arm number in the 
outer disc, and there is also confusion over whether there 
are 2 or 4 spiral arms ( Amaral Lepine|1997 Levine et al 


2006). The Be model has a weaker radial dependance, with 


the m = 2 mode being weaker in the inner disc compared to 
Bb. It does however have considerable power in the m = 5 
mode at around the Solar position. This is interesting as in 
many Galactic models the preferred structure is a 4-armed 
spiral, but with the addition of some small spur or minor Lo¬ 
cal arm which is always separate from the 4 main arms. The 
fact that the Be model is preferentially a 4-armed model 
throughout most of the inner disc, lightly 2-armed in the 
centre, and seemingly 5-armed near the Solar radius makes 
it an excellent candidate for replicating Galactic l-v features. 

The live disc-bulge model with the lightest disc (Bd) 
appears to be a clear departure from the Bb and Be models 
in that there is no clear dominant spiral mode visible by-eye. 
The structure of this disc appears to mimic a flocculent spi¬ 
ral galaxy such as NGG 4414 with multiple small scale arms 
and inter-arm features ( Braine et al.|19^ ). What structure 
there is favours a m = 2 mode in the inner disc, similar to 
the other models, and a m = 5 mode in the outer disc. 

The pitch angles have a weak correlation to disc mass, 
with an approximate value across all models of 21°. This 
value is within agreement of the values suggested by the 
literature the Milky Way (e.g. Levine et al. 2006), though 
somewhat on the larger end of most estimates. The pitch 
angle for the lightest disc was not well constrained due to 
the highly flocculent nature of the disc. In all cases the log¬ 
arithmic structure is best fit in the inner/mid-disc. 

The pattern speeds in Table cle arly decrease with ra¬ 
dius for all models (see also Eig. A4). The arms generated 
here are therefore material in nature, and are continuously 
sheared out and re-formed during rotation. The dominant 
arm number tends to be maintained, with arms being re¬ 
built from the remnants of the wound-up arm. The resulting 
pattern speeds range from around 40-35kms“^ kpc“^ in the 
inner disc to 30-20kms“^ kpc“^ in the outer disc, generally 
greater than those values suggested for the arms in previous 
work, centring around Dp 20kms“^kpc“^ (e.g. Dias & 


Lepine||2005 Gurhard||2011| Pettitt et al.||2Q14 ). 


In Eigurej^the evolution of Qs is shown from 0-1 Gyr for 
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Figure 3. Evolution of the Ba model, where the heaviest disc mass is investigated, showing the stellar (top) and gaseous 
(bottom). Large scale arm structures can be seen, yet appear rather irregular with many knee and kink features. 
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Figure 4. As Fig. but for the Bb configuration. Arms appear more smooth and regular compared to the Ba model. 
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ure 5. As Fig. [^but for the Be configuration. A significant number of small arm structures are now visible compared to Ba 
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Figure 6. As Fig. but for the Bd configuration. The arm structure is very small scale and fiocculent with no clear dominant spiral 
mode visible by-eye. 
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Model 

JV(i?in) 

^(I^mid) 

N{Rout) 

Dp(i^in) [kms ^ kpc 

Hp(i?mid) [kms ^kpe 

Dp(i7out) [kms ^ kpc 

«[°1 

Ba 

2 

5 

3 

38 

44 

33 

19 

Bb 

2 

4 

4 

34 

25 

21 

20 

Be 

4 

4 

4 

37 

27 

25 

23 

Bd 

2 

4 

5 

39 

35 

28 

22 

Db 

3 

3 

5 

N/A 

N/A 

N/A 

27 

Dc 

3 

5 

5 

N/A 

N/A 

N/A 

28 

Hb 

3 

4 

4 

42 

29 

33 

26 


Table 2. The number of arms, pattern speed and pitch angle for each model. The number of arms and pattern speed are shown are three 
different radii as they are both a function of radius. The three radii of calculation are = 5kpc, i^mid = 7.5kpc and i^out = lOkpc, 
corresponding to the inner, mid and outer disc. Boldfaced arm modes denote the dominant mode and the arm number corresponding to 
the a fit. 



Figure 7. Evolution of the Toomre Qs-parameter in the Be disc- 
bulge calculation as a function of radius initially and after 0.1, 
0.2, 0.3 and IGyr of evolution. 


the Be model. Initially Qs ~ 1 in the mid-disc, as dictated by 
the setup conditions, making the disc borderline stable. Over 
SOOMyrs Qs can be seen to slowly increase throughout the 
disc, raising to approximately 1.2 in the mid-disc region. At 
much later times, of the order of Gyr, Qs 2 implying the 
disc is highly stabilised. The spiral structure of the various 
discs tends to prevail up to a Gyr, but begins to smooth 
out and dissipate at later times. The longevity of the spiral 
arms is also seen to be directly tied to the resolution of the 


stellar disc (Fujii et al. 2011), with simulations with only 
1 X 10^ particles displaying clear spiral arms for less than a 
Gyr. The issue of spiral longevity is a standing problem, and 


not one we aim to investigate here (see instead Sellwood & 


Garlberg||1984| |Sellwood||2011| |Fujii et aL||2011| |Baba et ah 

2013||D’Onghia et al.||20]T~^ 


The rotation curves of each live disc-bulge model show a 
dispersion that increases with disc mass. Dispersion around 
the mean rotation curve ranges from di50kms“^ in Ba to 
zb20kms“^ in Bd. The dispersion in the rotation curve 
therefore appears directly related to the number of arms 
formed, which in unsurprising as the dispersion is also incor¬ 
porated into the value of m predicted by swing amplification 
{m oc 5]“^, and E determines the stellar velocity dispersion). 


3.2 Simulations and morphology of other models 


In addition to our fiducial live bulge and disc models we per¬ 
formed additional calculations to test how omitting a bulge, 
and adopting a live instead of fixed halo modify our results. 
This included two models with no bulge component and two 
different disc to halo mass ratios (Db and Dc), and a model 
with a live halo (Hb). In Figure]^ we show the Hb, Db and 
Dc models in stellar (top) and gaseous (bottom) components 
after 250Myrs of evolution. 

The structure of the Db and Dc models is similar to 
the Bb and Be in the outer disc, but irregular in the centre 
without the inclusion of a bulge, showing a strong m = 3 
mode in the inner disc. The dominant arm modes for Db 
and Dc are m = 3 and m = 5 respectively, though both are 
present in each model (Table [^. The pitch angle however 
gives very high values, with values of around 27° for the Db 
and Dc models, far outside values inferred for the Milky Way 
and at the high end of values seen in external galaxies. The 
pattern speed is difficult to determine, as the arm structures 
in the mid/inner disc wind up very quickly due to lack of 
support provided by the inner bulge. 


Due to the nature of the bulge-free system, the rota¬ 
tion curve for these models decays rapidly approaching the 
Galactic centre. This leads to the lack of an inner Lind- 
blad resonance, which is believed to encourage the growth 
of m = 2 modes (i.e. a bar) due to the lack of a Q-barrier in 
the inner disc (Binney & Tremaine 1987 Gombes et al. 1995 


Dobbs Baba|2Q14| . If we allowed the bulge-free models to 

evolve for the order of several Gyr then an inner bar compo¬ 
nent is seen to develop, but at the expense of arm structure. 


For the live halo-disc-bulge model, Hb, the morphol¬ 
ogy is not dissimilar to the static halo models, specifically 
appearing somewhere between Ba and Bb in terms of arm 
features. The difference in morphology between this and 
the Bb model is likely caused by the additional dynami¬ 
cal flexibility caused by the resolved halo, which may be 
under-resolved to perfectly match the analytical potential, 
as well as the added errors induced by the artificial radial 
truncation. There appears to be a dominant 3-armed pat¬ 
tern in the inner/mid disc, which dissipates in the outer 
disc where a m = 4 mode then dominates. The pattern 
speed for the m = 4 mode gives pattern speeds in the range 
30km s“^ kpc“^^ Qp ^ 45km s“^ kpc“^. 
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Figure 8. Structure of three separate Milky Way models in stellar (top) and gaseous (bottom) components after 250Myrs of evolution. 
The models are a live a disc-bulge-halo system (Hb, left) and a live stellar disc without an inner bulge with a moderate (Db, middle) 
and light mass disc (Dc, left). 


Model 

tbf [Myr] 

Robs [kpc] 

Lobs [kms 1] 

Fit stat. 

Ba 

197 

8.5 

215 

0.994 

Bb 

226 

8.5 

200 

0.857 

Be 

207 

7.0 

200 

0.833 

Bd 

207 

7.0 

205 

0.768 

Db 

235 

8.5 

200 

0.974 

Dc 

216 

8.0 

200 

0.931 

Hb 

216 

8.0 

205 

0.925 


Table 3. Results of fitting to lohs: Kibs? ^obs cind tbf for each 
model. A lower fit statistic signifies a better fit to the observa¬ 
tional data of Figure 


3.3 Simple kinematic l-v maps 

The process of fitting to /obs, K^bs, ^obs and tbf against the 
CO observational data of Figure was performed for each 
model, where a lower fit statistic signifies a better fit to 
the observational data of Figure The resulting best fit 
l-v maps are shown in Figure with parameters for each 
model given in Table (Dc is omitted from the figure as it 
differs little compared to Db). The bulge-free models (Db, 
Dc) have some of the poorest fit values, which appears to be 
due their lack of emission at high velocities in the inner disc. 
Some arm structures are seen in the outer disc, but these do 
not stray far from the local velocity (vios ~ Okms“^). 

The remaining models, all with a live bulge component. 


provide a variety of l-v features. The heaviest disc model 
(Ba) appears to have a large velocity dispersion in the inner 
disc. While it has an inner structure that is aligned similarly 
to observations, there is a large amount of material at high 
velocities and not matching the observed features. Coupled 
with the irregular arm structure seen in the top-down map, 
we conclude this (our heaviest disc) is a poor match to the 
Milky Way and do not produce any full radiative transfer l- 
V maps. Moderate-to-light mass discs (Bb, Be, Bd and Hb) 
provide a better agreement (and lower fit statistics). The 
lightest disc produces a near-uniform emission structure in 
the inner disc due to the flocculent nature of the arms. The 
Bb and Be models provide good representations of the Ca¬ 
rina and Local arms, while producing an Inner Ridge of the 
correct orientation. 


3.4 Best fitting models and radiative transfer l-v 
maps 

We created synthetic emission maps using the radiative 
transfer code TORUS for the models at the best-fitting time- 
frame found in Section 13.31 The fit statistic has been re¬ 
calculated from the TORUS models now that the intensity 
correctly takes into account emission and absorption effects. 
The fit statistic follows the same trends across the maps as 
when we used the simple map creation tool. The three best 
fitting maps (in order) are the Bd, Be and Bb models. In the 
case of the Be model there was a later time-frame that had 
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Figure 9. Best fitting l-v maps are shown from the simple fitting procedure to find lohs^ ^obs cind tbf for models Db, Ba, Bb, Be, 

Bd and Hb. Best fitting values are given in each panel, along with the value for the fit statistic used to constrain them. V^bs cind Rohs 
are in kms“^ and kpc respectively. The brightness of the emission has been scaled to the range seen in observations, whereas that in 
Figure [To] has been independently calculated. 


almost as good a fit (0.846 on the scale of the fit statistics in 
Figure]^, almost lOOMyrs later at 292Myrs, the maps for 
which we also include here. 


In Figurep^ we show our four best maps, from top to 
bottom: Bb, Bc(207Myrs), Bc(292Myrs) and Bd. In the left 
column we show the top-down gas distribution, and in the 
right the TORUS CO l-v emission maps, created at the best 
fitting values of /obs, I4bs and Rohs- In each frame we also 
label significant arm features, with the same nomenclature 
as used for the Milky Way (analogous to Fig. 25 in Paper 1). 


In the heavy disc (Bb, top panel) model the inner l- 
V features are smiler to the observed Inner Ridge, though 
not as clearly defined. The top-down map shows this is a 
combination of a far and near arm feature. The near arm 
(labelled Carina) appears weaker than the far inner arm 
(labelled SCC; Scutum-Centaurus-Crux). This allows for the 
reproduction of the Carina Arm in l-v space, which requires 
an arm to be very close to the observers location, while 
avoiding large amounts of emission at local velocities in the 
range |/| < 50° as Carina Arm passes the line-of-line to 
the Galactic Centre. This model also has a feature at very 
low velocities similar to the 3kpc-Expanding Arm, though is 
angled steeper in l-v space than observations. A Local Arm 
feature has been produced by a spur of the Carina Arm, 
lying very close to the red-cross in the top-down map. This 
has been suggested by other studies; that the Local Arm 
is in fact some spur or inter-arm structure, rather than a 


primary arm (though the picture is still not clear; Reid et al. 
Xu et al.|2013 ). The caveat of this model is the outer 


2009 


arm structure. Both the Perseus and Outer arms, while clear 
in the top-down map, are weak or incorrectly placed in l-v 
space. The Perseus arm appears at velocities too similar to 
the local values, making it nearly indistinguishable from the 


Local Arm in CO emission. The Outer Arm can barely be 
seen in emission, its presence only given away by a couple 
of dense pockets of gas. 

The second model (Be) has a slightly lighter disc, and 
appears a much better fit for many of the l-v features. The 
arms in the second and third quadrant are an especially good 
reproduction of observations. In Figure[^we show a zoom in 
of this region shown alongside the observational data. The 
Local, Perseus and Outer arm features are all reproduced 
and have comparable line-of-sight velocities. The emission is 
however still somewhat higher than that observed, a problem 
with all maps produced. The top-down map shows the Local 
and Perseus arms are a bifurcation of the same arm. The 
Carina hook structure is reproduced but is present at the 
incorrect longitude. The tangent point lies upon Vela (/ = 
100°) rather than the observed Carina tangent (/ = 80°). 
This feature can be made to match better by increasing 
Robs, but at the expense of the other arm features. The 
Inner Ridge of this model is somewhat poorer than the other 
models. There is a void of emission at approximately I = 20° 
vios = 80 km s“^ where clear emission is seen in observations. 
The incorrect reproduction of the Inner Ridge is due to the 
SCC Arm tracing a near-circular arc in the inner disc, which 
is seen as the steep straight line in the l-v map. In the other 
models the SCC arm clearly “winds”, i.e. has a non-circular 
shape, so is seen to be angled in l-v space. 

The second Be map, created at approximately lOOMyrs 
after the first, appears to be the best reproduction by-eye. 
The model shows the Carina, Perseus, Outer and Local 
Arms as well as an Inner Ridge that is aligned similarly 
to observations. This model offers the best reproduction of 
the Inner Ridge and Carina arm simultaneously. The Ca¬ 
rina Arm appears to branch away from the SCC Arm (the 
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Figure 10. Our four best-fitting CO radiative transfer l-v emission maps with their x-y counterparts. The models from top to bottom 
are Bb, Bc(207 Myrs), Bc(292 Myrs) and Bd. The top-down maps only show gas material that is seen in CO l-v space; that of the highest 
density. The cross indicates the observer’s position (which differs between models). SCC refers to the Scutum-Centaurus-Crux arm in 
the 4-armed paradigm of the Milky Way, also referred to in the main text as the Inner Ridge when viewed in l-v space. Arrows indicate 
locations of prominent features in l-v space. 
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Figure 11. Zoom in of the first and second quadrant arm features in CO from the 207Myr Be arm model (second row, Fig |10| ). 
Observational data is shown in the upper panel for comparison. As the synthetic has stronger emission on average we have increased the 
emission of the observational data slightly to be on the same scale as the synthetic map in this figure. 


source of the Inner Ridge) allowing it to be correctly placed 
in l-v space without causing spurious emission in the inner 
disc, a problem the symmetric fixed potential models persis¬ 
tently encountered. There also appears to be a 4-armed outer 
structure, with a strong 2-armed inner structure, which adds 
weight to the models in the literature that suggest a strong 
stellar 2-armed component with weaker 4-armed one in the 
gas/dust and young stars 
|2009| . Local material is 
Perseus Arm. The Perseus Arm itself is hard to differen¬ 
tiate from the Local and Outer Arm features, which is the 
main problem with the model. All arm structures in the sec¬ 
ond quadrant appear at too shallow velocities, implying Vohs 
is incorrect or that the model rotation curve is too shallow 
near the Solar Radius. 

Our final model has the lightest disc, and the lowest fit 
statistic. As with the previous model, there is a good repro¬ 
duction of the Inner Ridge, Carina Arm and Local material. 
The Perseus and Outer Arms appear too weak, and there is 
a significant amount of emission in the inner disc (|/| < 30°). 
While common to all models, this excess emission is espe¬ 
cially evident here due to the general flocculent nature of 
the spiral arms. There is no clear inner disc structure and 
the many smaller arm features in the inner disc are seen in 
molecular emission, appearing as a great swath rather than 
distinct arm features. The flocculent nature is also the cause 
of the weakness of the Perseus and Outer Arms beyond the 
Solar radius. It is surprising that the emission features can 
be well reproduced by a model with seemingly no clear dom¬ 
inant spiral mode (Fourier analysis indicates m ~ 5). While 
the fit statistic indicates this model is a good fit, this is likely 
because there is little emission seen in the incorrect place, 
coupled with the correct reproduction of the Inner Ridge, 
Carina and Local Arms. The arm features however seem 
too weak in the outer disc for this model to be a correct 
reproduction of our Galaxy. 

Overall, the fit statistic favours the live bulge-disc 
model with the lightest disc, Bd. However, the features ap¬ 
pear too flocculent in the outer disc, and so we favour the 


( Drimmel|2000| Churchwell et al. 


again formed by a spur off the 


second best fit model. Be. The Bd model provides an excel¬ 
lent fit to the strongest emission regions, that of the Inner 
Ridge, but the Be reproduces other arm features better. 


4 DISCUSSION 
4.1 General results 


In this paper we have shown that a live iV-body system rep¬ 
resenting the stellar component of the Galaxy can provide 
a good match to many of the observed molecular emission 
features. While not one single model produced all arms per¬ 
fectly, we believe that with enough initial seeds a match 
could be found. The Be model at the later time-stamp in 
particular compares to the Milky Way remarkably well, with 
a good reproduction of all arm features and inner emission 
structure. 

The values in Table [ 3 ] show that models with the lowest 
fit statistic (< 0.9) have a high arm number (3 < A" < 5). 
Models showing 2-armed structures were not readily pro¬ 
duced. In the case of the Ba model a 2-armed morphology 
was highly irregular and appeared to be buckling in the outer 
disc, beyond the Solar position. 2-armed models have been 
produced by studies in the literature, but tend to only be so 


when perturbed by some external body (Toomre & Toomre 
1972| Salo Laurikainen||2000| [Dobbs et al.||2010| [Struck 


et al. 


2011). The study of 


Martos et al. 


( 2004| ) proposed 
a 2-armed spiral potential that drives a 4-armed structure 
in the gas. However, in Paper 1 this model was tested and 
provided a poor match to the l-v data. 

The spiral arms produced here are transient in nature. 
Over the order of a Gyr the dominant arm number will be 
consistent for a given model, but the stellar and gaseous 
material for each arm reside in that arm until it is sheared 
apart by differential rotation of the disc. This is at odds 
with the spiral density wave theory that suggests arm fea¬ 
tures would prevail over long time periods, with material 
continually flowing into and out of the potential well of the 
arms. The fact that the transient arm structures presented 
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here are a better match than the steady waves of the po¬ 
tentials used in Paper 1 (see Sec. 4.2), and that standing 


spiral waves are yet to be reproduced in any numerical sim¬ 
ulation is evidence that the density wave theory may not be 
applicable to the Milky Way. 

There is however one key component absent from the 
models shown. The central molecular zone is very broad 
and featureless, a stark contrast to observations and pre¬ 
sumably a result of the lack of an inner bar. This is due to 
the Q-barrier caused by the inclusion of the bulge required 
to match the rotation curve of the Milky Way. Interestingly 
all models show a trailing m — 2 mode in the inner disc. It 
is possible this is the disc attempting to form a bar but is 
undermined by the increased stability inherent to the bulge 
dominated region. While we did manage to produce bars in 
some of our preliminary simulations, they required a long 
term evolution of the stellar disc of the order of many Gyrs 
by which point the strength of the arms is greatly reduced. 
Additionally the lack of a bulge means the inner rotation 
curve is a poor match for the Milky Way. The generation of 
Wbody bars in Milky Way like models with a co-existing 
strong spiral structure will be the subject of a future study. 


4.2 Comparison to calculations with fixed stellar 
potentials 

In Paper 1 we performed calculations using a number of 
fixed analytic potentials to represent the arm and bar com¬ 
ponents of the Milky Way. A 4-armed model was favoured 
to match all the arm features, with pitch angles around 15°, 
pattern speeds of 20km s“^ kpc“^ and 50-60kms“^ kpc“^ for 
the arms and bars respectively. While many arm structures 
could be reproduced, the regular structure made it impossi¬ 
ble to create a perfect match. For instance, the Carina arm 
was impossible to correctly reproduce without placing very 
strong local emission in the inner Galaxy. 

It appears that using a live-stellar distribution provides 
a much better match for Galactic l-v structure than the 
fixed analytic potentials in many regards, clearly seen by 
comparing Figure 25 of Paper 1 and Figure here. The ir¬ 
regular arm structures created by the live stellar system are 
able to match emission features simultaneously, such as the 
Carina arm and Inner Ridge, where symmetric logarithmic 
spiral arms could not. For example, in the models of Pa¬ 
per 1 the second quadrant could be fit by moderate to large 
pitch angles, whereas a much smaller value was needed to 
fit the Carina feature. The arm numbers of the best fitting 
models are similar to those suggested previously, favouring a 
4-armed gas structure to best match the observations. The 
pitch angles are somewhat higher than found in Paper 1 


the case that the Galaxy has no fixed pattern speed, with 
spiral arms that are also material in nature. Some external 
galaxies have also been observed to have a pattern speed 



without the arms dissipating. Pattern speeds appear to be 
a function of radius in nearly all cases, with mean values for 
arms ranging from 20km s“^ kpc“^< Qp < 40km s“^ kpc“ 


also similar to values found in other studies ( Grand et al. 
2012 Baba et al.||2013 ) but higher than those found when 
using analytic potentials. The arms appear material, unlike 
the steady density waves implied by theories, and material 
tends to reside in the arms until they shear apart. It may be 


that decreases with increasing radius (|Meidt et al. 
Speights Westpfahl|2012 ). 


2008 


The total strength of CO emission is more in line with 
that seen in observations (i.e. weaker) than the models of 
Paper 1, which tended to create extremely bright arm emis¬ 
sion structures throughout the arms. The increased amount 
of inter-arm structures allows for lower strength emission 
features, whereas analytic potentials tended to create one 
strong swath of emission tracing the potential minimum. 
The fit statistic for the emission maps of models is calcu¬ 
lated in exactly the same way as in Paper 1 and the best¬ 
fitting models here provide a systematically better fit than 
those with fixed potentials (best fitting values give ~ 1.05 
for potentials and 0.95 for live discs). 

There are some drawbacks to this method however. 
Each model has far too much emission in the inner disc. 
The fixed potential calculations effectively had a hole in the 
inner disc, which resulted in a large dearth of emission at 
high velocities inside of |/| < 20°. High gas density was still 
seen in the inner disc, but was solely aligned on the X 2 orbits 
of the bar. As such we believe either an inner bar structure 
is needed to sweep up molecular material in the inner disc, 
or that gas density is greatly reduced by some other mech¬ 
anism. 


5 CONCLUSIONS 

In this study we have shown simulations of the stellar and 
gaseous components of the Milky Way. Different spiral mor¬ 
phologies were formed in the stars and gas, with arms ap¬ 
pearing transient and material in nature rather than as den¬ 
sity waves. The arm number is seen to increase as the disc 
to halo mass ratio decreases, with arm numbers found to 
range from 2 ^ iV ^ 5 similar to those predicted from 
swing amplification theory. We perform fits to logarithmic 
spiral features, finding pitch angles in of 18° < a < 25° 
and pattern speeds in the range of 20kms“^ kpc“^< Dp < 
40kms“^kpc“^ which decreases with radius rather than 
maintaining a constant value. Both pattern speed and pitch 
angle are within the range of values inferred from observa¬ 
tions, though in the higher part of this range. 

Using the molecular gas in these simulations we then 
created synthetic l-v emission maps. A simple method is 
used to find a best-fitting time-frame and observer coor¬ 
dinates, which are used to reject some outlying models 
and provide input parameters for the full radiative trans¬ 
fer maps. We find moderate mass discs (model Be; Md = 
3.2 X 10^°M©) with a live bulge-disc component provide a 
very good match to the observations, with 4-armed spiral 
patterns that reproduce many of the arm features. These 
arms provide a better fit than those using fixed potentials 
and provide a lower fit statistic. The arm features of the 
Milky Way are thus found to be best-fit by a dynamic and 
transient disc, displaying a predominantly 4-armed pattern 
in the gas with a pitch angle of approximately 20°. Some ob¬ 
servational studies propose a clearer 2-armed structure (e.g. 
Weaver|1970 Driminel|2000 Steiman-Gameron et al.|2010 ), 
whilst some are more favoured towards a 4-armed structure 
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(e.g. 

Georgelin V Georgelin|1976| Russeil|2003 

Levine et al. 

|2006 

), though we note again that our models tend to display 


2 arms in the centre of the disc. 


The main features of the Milky Way appear to be re¬ 
producible with relatively simple physics. By only taking 
into account the stellar gravitational held we can create a 
strikingly similar match to the observed features. This im¬ 
plies that, for most of the Milky Way disc at least, the effect 
of feedback, self-gravity, the bar, and perturber interactions 
can be considered minimal compared to the stellar gravity. 

A clear improvement to the work in this and Paper 1 
is the inclusion of additional physics. Aside from hydrody¬ 
namics and chemistry, our calculations are relatively simple. 
The logical next step would be to take the best-htting mod¬ 
els from the analytic potentials and the live stellar systems 
and include gas self-gravity, stellar feedback or even Galac¬ 
tic scale magnetic fields. It would also be interesting to see 
whether tidal forces induced by an object of the order of the 
mass of the LMC could induce a spiral structure resembling 
what is seen in l-v Galactic observations. 
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APPENDIX A: QUANTIFYING ARM 
MORPHOLOGY 



6 [rad] 


Figure Al. The stellar component of the Be model binned into 
R — 6 space and normalised by azimuthally averaged surface den¬ 
sity. The diagonal features show location of spiral arms that are 
detected by the Fourier decomposition. 


pervious studies ( Dobbs et al.|2010 Eujii et al.|2011 Grand 


et al.|[MT^ D’Onghia et al.||2013 ). An important caveat is 

that the arm shape is assumed logarithmic and periodic, i.e. 
with each arm separated by 2ii/N. The particles are first 
binned into a R-0 grid which is then normalised by the sur¬ 
face density of the disc to ensure that inner arm structures 
do not dominate the fit. This gives the surface density con¬ 
trast at a given radius: r|i^(^) = X(R, 0)/T^{R). An example 
of which for the Be model is shown in Eigure fAl] 

A Eourier transform is then performed on r|i?(^) at dif¬ 
ferent radii, the results of which for the stellar and gaseous 
components of our four main models (using a live disc and 
bulge) are shown in Eigure [A^ The dominant modes clearly 
appear to increase with decreasing disc mass, with the m = 2 
mode dominating the Ba model, and the m = 5 mode gain¬ 
ing significant power in the Be and Bd models while appear¬ 
ing negligible in Ba and Bb. The strongest Eourier mode is 
taken to be the arm number, though in some instances there 
are two conflicting modes due to the variation as a function 
of radius. We take the m = 4 mode in the Bb model, as it 
is stronger in the stellar component and near the solar po¬ 
sition. Eor Be the m = 4 mode is also selected, though the 
m = 5 mode displays significant power over a wide radial 
range in power components. Eor the Bd model we select the 
m = 5 mode. The m = 2 mode is also strong, but only in the 
inner disc, whereas the m = 5 mode is strong in the outer 
disc yet still peaks in the inner disc. 

We then extract the R-0 position of the dominant mode. 
Depending on the morphology, this may only be done across 
a certain radial extent, such as the range in which that mode 
dominates (e.g. from 2kpc ^ R ^ 6kpc for the m — 2 mode 
in the Bb). The pitch angle, a, of a logarithmic spiral arm 
is linked to the values of the constant B in the equation 


0 = /(R,a) = RlnR + C 


(Al) 


The arm number (iV), pitch angle (a) and pattern speed 
(Qp) are determined using methods similar to those used in 


where the pitch angle is a = arctanR and C defines the 
azimuthal position of the arms. This equation is then fit to 
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the relevant mode using a downhill Nelder-Mead simplex 
algorithm from the SCIPY PYTHON package and minimising 
a chi-squared like statistic. 

Figure [A3] shows the fits to the dominant stellar compo¬ 
nents of the models in Figure [A^ Black points trace out the 
arm and the best fitting logarithmic spiral arm is shown by 
the red line. The logarithmic spiral appears a very good fit 
to the arm features in our models, though the weak nature 
of the Bd model arms leaves a larger margin of uncertainty. 

We then take subsequent time-steps and perform a sim¬ 
ilar analysis, this time fitting to only the dominant spiral 
mode determined from the main time-stamp. Then by sim¬ 
ply calculating the offset between arms at different epochs 
the pattern speed can be calculated as a function of ra¬ 
dius. In Figure [A4| we show the pattern speed measured for 
the models shown in Figures |A2| and |A3| The individual 
points are calculated over a range of ?^40Myrs, with error 
bars showing the maxima and minima at lOMyr intervals 
within this window. In all instances a clear decrease of pat¬ 
tern speed with radius can be seen, indicating the arms act 
as material, rather than wave-like features. 

This paper has been typeset from a T^]X/ BTgX file prepared 
by the author. 
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Figure A2. The power of individual spiral modes in the range 1 ^ m ^ 6 for the gaseous and stellar components of models Ba, Bb, Be 
and Bd. Similar modes dominate in the gas and stars with higher arm modes dominating the lower mass discs. 
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Figure A3. Best fitting logarithmic spiral model (red line) to the dominant modes given by |A2| for our live disc-bulge models (black 
points). The best-fitting pitch angles are given in the top left of each panel. 
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Figure A4. Pattern speeds as a function of radius for the arms fit in Figure [A^ Each point is calculated from a time interval of 40Myrs, 
with error bars showing the maxima/minima in pattern speed measured at lOMyr intervals in the 40Myr window. 

























